Magnetic Properties of a Bose-Einstein Condensate 
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Three hyperfine states of Bose-condensed sodium atoms, recently optically trapped, can be 
described as a spin-1 Bose gas. We study the behaviour of this system in a magnetic field, and 
construct the phase diagram, where the temperature of the Bose condensation Tbec increases 
with magnetic field. In particular the system is ferromagnetic below Tbec and the magnetization 
is proportional to the condensate fraction in a vanishing magnetic field. Second derivatives of 
the magnetisation with regard to temperature or magnetic field are discontinuous along the phase 
boundary. 



Recently experimentalists confined three / = 1 hyper- 
fine states of sodium in an optical trap W. This was 
impossible to do in a magnetic trap, which confines only 
one hyperfine state [^j. One can ask what are the prop- 
erties of an / = 1 Bose-Einstein condensate in a mag- 
netic field. In a weak magnetic field H the electron s 
and nuclear i magnetic moments precess around f and f 
precesses around H. A simple calculation |^ gives that 
the average of the projection of the electron spin on the 
direction of the magnetic field (z-direction) is: 

[f{f + l)-^{^ + l)+J{3 + l)]■ (1) 



2/(/ + l) 



Sodium has a nuclear spin « = | and electron spin s = ^ , 
and the electron is in an Z = state, so that j = s. 
Substituting these values in Eq. (pT), we obtain: 



1 



(2) 



The magnetic moment of the nucleus is negligible in com- 
parison to that of the electron. Therefore, the system is 
equivalent to a spin-1 boson, with magnetic moment 4 
times less than that of the electron. The projection of 
the magnetic moment on z-direction is: 



o l^B 



(3) 



where /is is the Bohr magneton. As the hyperfine inter- 
actions are of the order of GHz, magnetic field much less 
then GHz/^B wlOOGs is small. 

As the densities of the Bose condensates of atomic 
vapours are very small, one can treat them as an ideal 
gas. We shall return to the validity of this assumption 
below. 

An immediate observation is that the system is ferro- 
magnetic 1^, when Bose-condensed, as the condensate is 
aligned with the external magnetic field however small 
it is. Another observation is that the external magnetic 
field increases the temperature of the Bose condensation. 
A spin-s Bose gas has a critical temperature which is 
(2s -I- 1)^/'^ times less than a spin-0 one at the same den- 
sity, due to the 2s -I- 1 spin states available for excitations. 
At high magnetic fields the gas is completely spin polar- 
ized and therefore equivalent to a spin-0 Bose gas, and 



has the same critical temperature as the latter. The aim 
of this paper is to construct the phase diagram of the sys- 
tem, to find the condensate density and magnetization as 
a function of temperature and magnetic field and investi- 
gate singularities of the latter along the phase boundary. 

The number of particles in excited states is given by a 
standard formula |||] , modified to include 2 extra spin 
states: 



imatT)i 
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(4) 



Here 



il^B/2)H 
T ' 



(5) 



fi is the chemical potential and mat is the mass of the 
bosonic atom. The integrals in Eq. can be rewritten 
following Q, Q as 



with F defined as: 



T 2^ 



-Fdx), 



(6) 



(7) 



n=l 



which reduces to the Riemann function for x 
0: Fa{0) = Cl"^)- Let us introduce new variables: 



To - Tbec{H - 0) 



27r 



(3C(i))2/3 mat 



2/3 
"^0 ' 



t = T/To, b- 



Tn ' 



(8) 



where no is the total particle density. 
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FIG. 1. Contour plot of the Bose-condensate density in 
/ = 1 sodium atoms (spin-1 bosons) in a temperature 
[t = T/To) - magnetic field (6 = i||^)plane. The solid (and 
uppermost) line is the boundary of the phase diagram. The 
dotted lines are contours for n = UBEc/no = 0.1,0.2, . . . ,0.9 
from top to bottom. Contour lines pile up near & = due to 
the square root singularity in the condensate density (see the 
text) . 

In the notation of Eq. (|) and Eq. (|), Eq. (|) can be 
rewritten as 

^ejno = (Fs (a) + F3{a + y) + Fs {a + 2y)) . 

(9) 

Here, as in the familiar spin-0 boson case a = 

when nex{a = 0,t,y)/nQ < 1 and a is determined from 
nexio,,t,y)/nQ — 1 othervise. 

The temperature, tc of BEC is obtained by equating 
the number of particles in excited states Uex to the total 
number of particles: n^x = n-o ^t zero chemical potential. 
Substituting this in Eq. (^) we obtain: 
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i = + + (10) 

This gives the phase diagram in a parametric form. In 
order to obtain it explicitly, one computes numerically 
tc{y) from Eq. (|l^ ) and then finds bc{y) = y/tc{y). To 
obtain the condensate density, n — ubec /i^a, contour 
lines, we replace 1 in the left hand side of Eq. HO with 
1 — n. Then for given y we get f„ = tc{l — nfi^ and 
bn = y/tn{y) = i'c{y)/il — n)-^^^. Therefore the conden- 
sate density contour lines are obtained from the phase 
transition line by a rescaling. The results are shown in 
Fig. ^ For small or large b one can find the asymptotic 
behavior of Eq. (p^): 

i^^H.l(i±^^V6,when6«l, 



tc ^ 32/3 - ^Y7^ exp(-V3'/'), when b » 1. (11) 

When 6 CX3 all atoms are in the mf — I state which 
is equivalent to having a spin-0 Bose gas at the same 
density, resulting in an increase of tc by a factor of 3^^^. 

The condensate fraction n = nBEcIno = 1 — Uex/nQ, 
computed using Eq. (||) as a function of magnetic field 
or temperature is shown in Fig. ^. Asymptotically, the 
condensate density is n{b = 0) = l — t^^'^ and n{b) — n{b — 
0) t\/b for i < 1 and 6 < t. For t > 1, n is linear in 
b-bc. 

To calculate the magnetization per particle, m, we note 
that (in units of ^): 

TO — {n{mf = 1) — n(mf = —l))/nQ. (12) 

This expression is inconvenient because there are both 
condensate and excited atoms in the mf = 1 state. Tak- 
ing into account that n{mf = 1) -I- n{mf = 0) -t- n{mf — 
— 1) = no, one has to = 1 — n{mf = 0) — 2n{nif = — 1), 
or: 




0.2 0.4 0.6 0.8 1 



b 

FIG. 2. Fraction of condensate, n = UBEc/no, as 
a function of temperature {t = T/Tq) and magnetic 
field {b = -^i^^p^). (a) n{t) for several values of 
b — 0,0.1, 1, 10, from left to right, (b) n(6) for several val- 
ues of t = 0.5, 0.8, 1, 1.2.1.4, from left to right. 
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FIG. 3. (a)Magnetization, m , as a function of temper- 
ature {t = T/Tq) for several values of the magnetic field 
(6 = LubME. ^ 0,0.01,0.1,0.3,1,3, from left to right). 
Dashed line is m{bc{t),t). (b) derivatives of the m{t) curves 
in (a) with regard to temperature. The cusp is at the critical 
temperature. 



Asymptotically: 



m = 1 — t2 



(14) 



when b = 0+, that is the magnetisation is proportional 
to the condensate density in a vanishing field, and 



Vb, 



(15) 



when i = 1, i.e. at the phase transition. 

When b is finite, the magnetization is no longer zero in 
the normal phase, however BEC appears in a cusp in the 
derivatives (^){, and (^)j, see Figures || and ^, where 
t hes e curv es w ere computed numerically, using equations 
(Al) and (|A4|) . In other words, second derivatives of the 
magnetization are d isco ntinuo us a long the phase bound- 
ary (see Equations (A3) and (A6)). 

We have studied the magnetization of the Bose- 
ferromagnet in a vanishing but non zero magnetic field. 
A question remains what happens for strictly zero field. 
In the case of Bose condensate of N spin-1 particles pro- 
jection of the magnetization per particle, m (in units of 
particle's magnetic moment) on the z-axis is equal to 
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FIG. 4. (a)Magnetization, m , as a function of mag- 
netic field (6 = liLsipil) for severall values of temperature 
{t = T/To = 0.8, 1, LI, 1.2) from left to right). Dashed line 
is m{b,tc{b)). (b) derivatives of the m{b) curves in (a) with 
regard to magnetic field. The cusp is at the critical magnetic 
field. 



TO, = (A^(l) - N{-l))/N. 



(16) 



Ah values of iV(l), iV(0), A^(-l), which satisfy the obvi- 
ous restriction 



N{l) + N{0) + N{-1) = N, 



(17) 



have equal probability, = 1 is achieved only when 
iV(l) = N,N{0) = 0,7V(-1) = 0, while = can 
be achieved by many combinations: A^(l) = 0,N{0) — 
N, N{~1) = 0, or N{1) = 1, iV(0) = iV - 2, iV(-l) = 1 
and so on. One can show by a generalization of this ar- 
gument that the probability distribution is 



P,{m,) = 1- 



(18) 



To find how the absolute value of to is distributed we 
notice that the magnetization is macroscopic and there- 
fore can be treated classically. Now if we know the dis- 
tribution of its projection on the z-axis Pz{mz), we can 
find the distribution of the absolute value P{m). Tak- 
ing into account that the projection of a vector of unit 



3 



length uniformly distributed on a sphere is uniformly dis- 
tributed between —1 and +1, we can write: 

Pzi^nz) — I dniP{m)- j dx6{mx — rriz) = 



2 7_i 



dmP{m)/{2m). (19) 



By differentiating this equation with respect to we 
get: 



P{m) = -2mP^(m), 



or, taking into account Eq.(18) 
P{m) — 2m 



(20) 



(21) 



These distributions are shown in Fig. |5|. Let us compare 
this with classical systems. For instanse in the familiar 
case of an Ising ferromagnet below the phase transition 
temperature the probability distribution of the magneti- 
zation consists of two one-half delta functions at positive 
and negative values of the magnetization. Applying a 
magnetic field changes this distribution to a single delta 
function. For the case of the Heisenberg ferromagnet 
magnetization can point in any direction, however the 
probability distribution of the absolute value of the mag- 
netization is a delta function, an infinitely small externall 
magnetic field will select the direction of the magnetiza- 
tion. In the case of a Bose-ferromagnet, an infinitely 
small applied magnetic field will not only select the di- 
rection, but also change the form of the probability dis- 
tribution fimction of the magnetization from that of Fig. 
||(b) to delta-function. 

The experiments are done in a trap, which can be mod- 
eled by a weak external potential. It was shown that 
interatomic interactions, however small, have a strong 
effect on the properties of the condensate, when it is in 
an external potential. The effect is that the Bose con- 
densate is not in the single-particle ground state of the 
potential, which has small spatial extent, but its wave 
function is spread to balance the interatomic repulsions 
and confining potential. The density of the condensate 
is then determined from this balance. Once it is deter- 
mined, it is a reasonable approximation to treat the gas 
as ideal at that density, as if it was in a square box, which 
is in the optical trap ~ lO'^^/cm'^. This makes < = 1 in 
our figures to be of the order of /iK and 6 = 1 of the 
order of O.lGs (which is also much less than the weak 
field requirement of lOOGs mentioned in the begining of 
the paper). At this density spin-spin interactions are a 
few orders of magnitude less than the spin energy in a 
magnetic field, and can therefore be neglected. A disad- 
vantage of this low density is that the magnetization is 
very small and difficult to observe directly. However one 
can meassure the populations of the hyperfine states |^] 
and calculate the magnetisation from Eq. (O). 
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FIG. 5. (a) Probability distribution of the projection on 
arbitrary axis of the spin-1 Bose condensate's magnetization 
in a zero magnetic field. Magnetization is in units of the mag- 
netization of completely spin-polarized condensate, (b) Same 
as (a), but for the absolute value of magnetization vector. 



Exchange interactions have been found to be impor- 
tant j|] in a related problem of a charged spin-0 Bose 
gas |10 (which is diamagnetic when Bose-condensed) . In 
our case they also have a profound effect and change the 
order of the transition from second to first. The previ- 
ous discussion was of ideal Bose gas in an external mag- 
netic field. As we argued, fields of the order of 0.1 Gs 
(which are also of the order of the magnetic field of the 
Earth) will dominate over exchange interactions, so the 
effect of the latter may not be observed. However, if 
the external field is screened down to 10~''Gs they may 
play a role. We shall now address the case of a Bose 
gas with ferromagnetic exchange interactions. Without 
ferromagnetic exchange the system undergoes a conven- 
tional Bose condensation, while without Bose statistics it 
undergoes a conventional ferromagnetic transition. This 
system is characterized by two parameters: the Bose con- 
densation temperature The, when ferromagnetic exchange 
is neglected, and the ferromagnetic transition tempera- 
ture Tf, when the effects of the Bose statistics are ne- 
glected. The behaviour of the system will depend on the 
ratio K = Tf /T^e ■ The details and the complete phase 
diagram will be reported in a future publications, here 
we shall give an essential result. If k > Kc, where Kc 
is close to 1, that is in the case of strong ferromagnetic 
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interactions, upon cooling, the system undergoes first a 
ferromagnetic transition and then at a lower tempera- 
ture a Bose condensation. Both transitions are of second 
order. If k < Kc then two transitions merge into a sin- 
gle Bose-ferromagnetic phase transition, which is of first 
order. Estimates show that for the condensates, experi- 
mentalists are currently working with, k << Kc- 

A few final remarks. Since in a vanishing field the 
condensate carries magnetization while the normal com- 
ponent does not, superflow is accompanied by a trans- 
fer of the magnetization. In particular, second sound 
should be accompanied by magnetization waves. The 
Einstein-deHaas effect should also be considered. When 
before condensation atoms are in equal proportions in 
•nif = —1,0,1 states, condensation will turn them all 
into an TO^ = 1 state, which should result in the creation 
of a single vortex. 

A spin-1 Bose gas is ferromagnetic when Bose con- 
densed. The reason for this is the Bose-Einstein statis- 
tics, and not the exchange interactions, like in conven- 
tional ferromagnets. Similarly, a charged Bose gas would 
behave like a superconductor JTot . 

The authors acknowlege helpful discussions with J.T. 
Liu, H.C. Ren and J. Stenger. This work was supported 
by the Department of Energy Contract No. DE-F602-88- 
ER13847. 



APPENDIX: 



From Eq.h3 one obtains: 



dm 
'dt 



t2 



-F.{a + y) + Fi{a + 2yy 



|(Fi(a + y) + 4i^i(a + 2y))- 

1 / OcL 

-iF,ia + y) + 2F,ia + 2y))ti^-£ 
Where (ff)j, can be found from Eq.(^ to be: 



(Al) 



da 

m 



-{F.{a) + Fs^{a + y) + F^{a + 2y))+ 
y{F^{a + y) + 2Fi{a + 2y)) 
Fi{a) + Fiia + y) + F.ia + 2y) 



(A2) 



Now (^)^ can be foimd numerically from Eqs. (Al), 
(|A2|), with a from Eq.M. To find the discontinuity 
in the at the phase transition we note that when 



0*2 

'da' 



a = 0: (If)j^ = (Cf.Eq.(|A2|), taking into account that 
Fl(0) — oo). Therefore the discontinuity comes only 
from the second derivative of a. We find for the jump: 
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d^m Fi{y) + 2Fi{2y) 



6<(i)Vi 
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F^i2y)) + yiF^iy) + 2F^{2y)) 



(A3) 



p2 1 

This has assymptotics ~ b^2 for 6 <C 1 and 

A^ ~ exp(-V3^) for b > 1. 

Derivatives with regard to the magnetic field can be 
done analogously, leading to: 



dm 
'db 



t2 



3C(i) 



F^ia + y) + AF^{a + 2y)y 



da 



da 
db 



iF.ia + y) + 2Fiia + 2y))t[ — 



Fiia + y) + 2Fi{a + 2y) 
'FUa)+FUa + y)+FUa + 2yy 



A- 



F^iy) + 2F^{2y) 



db^ 



6K(|)Vi 



This has assymptotics A 



d m 



and A|^ 



-(3t -tf for 35 



t < 1. 



(A4) 



(A5) 



(A6) 



for i - 1 < 1 



o2 

One can show similarly that is also discontinuous 
along the phase boundary. 
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